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We present a class of generally covariant nonlocal gravity models which have a flat-space general 
relativistic (GR) limit and also possess a stable de Sitter (dS) or Anti-de Sitter (AdS) background 
with an arbitrary value of its cosmological constant. The nonlocal action of the theory is formulated 
in the Euclidean signature spacetime and is understood as an approximation to the quantum effec- 
tive action (generating functional of one-particle irreducible diagrams) originating from fundamental 
quantum gravity theory. Using the known relation between the Schwinger-Keldysh technique for 
quantum expectation values and the Euclidean quantum field theory we derive from this action the 
causal effective equations of motion for mean value of the metric field in the physical Lorentzian- 
signature spacetime. Thus we show that the (A)dS background of the theory carries as free propa- 
gating modes massless gravitons having two polarizations identical to those of the Einstein theory 
with a cosmological term. The on-shell action of the theory is vanishing both for the flat-space and 
(A)dS backgrounds which play the role of stable vacua underlying respectively the ultraviolet and 
infrared phases of the theory. We also obtain linearized gravitational potentials of compact matter 
sources and show that in the infrared (A)dS phase their effective gravitational coupling Goff can 
be essentially different from the Newton gravitational constant Gjv of the short-distance GR phase. 
When Geff S> Gn the (A)dS phase can be regarded as a strongly coupled infrared modification 
of Einstein theory not only describing the dark energy mechanism of cosmic acceleration but also 
simulating the dark matter phenomenon by enhanced gravitational attraction at long distances. 

PACS numbers: 04.60. Gw, 04.62. +v, 98. 80. Bp, 98.80.Qc 
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I. INTRODUCTION 

Central problem in the attempts to build a modified 
gravity theory as a model of observable cosmic accelera- 
tion [ij] consists in elimination of ghost instabilities that 
usually make any model physically inconsistent. Com- 
bined with the necessity to retain the general relativistic 
limit, this makes the simplest appropriate version of gen- 
eral relativity (GR) - explicit cosmological term - very 
robust against possible attempts to modify it, either by 
introducing higher-derivative terms or relaxing its diffeo- 
morphism symmetry. 

On the other hand, there is fine tuning problem associ- 
ated with the hierarchy of the fundamental Planck scale 
vs the cosmic acceleration scale and also the cosmic coin- 
cidence problem - order of magnitude equality of dark en- 
ergy (DE) and matter contributions (including dark mat- 
ter (DM)). These problems serve as a strong motivation 
to go beyond introduction of explicit cosmological term 
and give rise to numerous models (like R + B? / K mod- 
els [2|, quintessence [3|1, /(i?) models [4 L non- minimally 
coupled matter fields [5| , brane theories 6[ , massive grav- 
ity 0,[1|, nonlocal cosmology etc.). However, in 
this or that way fine tuning is creeping into almost all 
of these models. Modulo certain exceptions [5], most of 
them in fact look as a sophisticated way to incorporate 
into their action in addition to the Planck scale the hori- 
zon scale (whether it is explicit cosmological constant, 
graviton mass of massive gravity [7[, multi-dimensional 
Planck mass in braneworld theories or the crossover scale 
in brane induced gravity models jg], etc.). 

To circumvent this difficulty one could adopt another, 
perhaps more promising, line of reasoning. If a con- 



crete fixed scale incorporated in the model is not sat- 
isfactory, then one could look for a model that admits 
cosmic acceleration scenario with an arbitrary scale. Its 
concrete value compatible with observations should arise 
dynamically by the analogue of symmetry breaking to 
be considered separately. Even this very unassuming 
approach is full of difficulties, because modified gravity 
models featuring this property (like unimodular gravity 
[13], /(i?)-gravity, etc.) generally violate some of its con- 
ventional symmetries, have additional degrees of freedom 
and might lead to ghost instabilities. Thus we get back 
to the central problem in the modification of Einstein 
theory, mentioned above, which is especially very actual 
under the requirement to preserve general covariance of 
the model. 

Here we present a nonlocal infrared modification of 
Einstein gravity theory briefly reported in [l^, which is 
likely to implement the above approach. It will be based 
on the realization of the old idea of a scale-dependent 
gravitational coupling - nonlocal Newton constant [T^ 
Il7| - and will amount to the construction of the class 
of diffeomorphism invariant models. These models are 
compatible with the general relativistic (GR) limit and 
generate a stable de Sitter (dS) or anti-de Sitter (AdS) 
background with an arbitrary value of its effective cos- 
mological constant A. Repeating the above motivation 
again, arbitrariness of A will be a manifestation of the 
fact that, to resolve such issues of DE as fine tuning and 
cosmic coincidence, the scale of A cannot be encoded in 
the action of the theory, but rather should arise dynam- 
ically by the analogue of symmetry breaking. 

In addition to fine-tuning argumentation of the above 
type, the driving force of our approach will also be the 
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aesthetical motivation to consider as a source of DE a 
purely metric sector of the theory. No special matter 
fields like quintessence [1], dilaton ^ or Lorentz break- 
ing khronon Ts] will be assumed to exist. In distinction 
from local f{R) models this will be a nonlocal generally 
covariant modification of the metric sector of Einstein 
theory. Also, it will go beyond the f{R/0) models of 
|9l4ll| by involving a very nonlinear and nonlocal depen- 
dence on all components of Ricci and Riemann tensors. 
Finally, though our model stems from the idea of a non- 
local gravitational coupling "constant" motivated by 
the idea of degravitation of the vacuum energy (l7j . no 
degravitation of matter will be considered, because the 
source of the effective cosmological constant will be en- 
tirely the metric sector of the theory. 

Although the main goal of this approach was the 
achievement of DE phenomenon, the serendipitous na- 
ture of the resulting model will appear as an unexpected 
bonus - dark matter simulation also generated by the 
metric sector of the theory. Whatever speculative is the 
interpretation of this effect, it arises as an enhancement 
of the gravitational attraction of the ordinary (non-dark) 
matter at large distances, which is quite opposite to the 
degravitation mechanism. 

The plan of the paper is as follows. In Sect. II we formu- 
late the four-parameter family of nonlocal gravitational 
models and give a brief summary of their properties - 
existence of a stable flat-space background correspond- 
ing to their GR limit and an alternative ghost- free (A)dS 
background with an arbitrary value of the effective A, 
corresponding to their infrared DE phase. Here we also 
give the expression for the linearized gravitational poten- 
tial of matter sources in this (A)dS phase, which has a 
DM interpretation. In Sect. Ill within a flat-space pertur- 
bation theory setup we give a motivation for the nonlocal 
structure of their action, based on the generally covari- 
ant concept of the scale-dependent gravitational coupling 
and the requirement of stability with respect to ghost 
modes. In Sect. IV we discuss treatment of nonlocality 
under the assumption that the suggested nonlocal action 
is in fact a certain approximation for the quantum ef- 
fective action obtained from some fundamental quantum 
gravity theory. Here we employ the relation between the 
Schwinger-Keldysh technique for expectation values [l^ 
and the Euclidean quantum field theory to show how a 
nonlocal effective action leads to causal retarded nature 
of nonlocal equations of motion for mean fields |2lj . We 
begin with this relation within the flat-space perturba- 
tion theory setup, but then show that in the cosmological 
context this setup fails unless we modify nonlocal struc- 
tures of the effective action by extra non-polynomial de- 
pendence on the spacetime curvature. In Sect.V we es- 
tablish this dependence in the action (thus recovering its 
final nonlocal form introduced in Sect. II) and prove the 
existence of the (A)dS background with an arbitrary A. 
Stability requirement for this background requires knowl- 
edge of the quadratic part of the action. This represents 
a technically simple but very lengthy calculational chal- 



lenge accomplished in this section and paper appendix. 
In Sect. VI we show that this (A)dS background carries 
as free propagating modes a massless graviton with two 
polarizations and also obtain the linearized gravitational 
potentials of matter sources, which is likely to simulate 
the effect of DM. Finally, in Sect. VII we discuss GR and 
(A)dS phases of our theory and the extent to which it 
might match with the cosmological concordance model 
and also its prospective nature in other ramifications of 
quantum gravity theory. 



II. SUMMARY OF RESULTS 

The action we will be interested in reads in its full 
complexity as the following nonlocal diffeomorphism in- 
variant functional of the spacetime metric g^i/,^ 
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Here G^i/ — R^^ — \g^yR is the Einstein tensor, the hat 
denotes the matrix acting on symmetric tensors and we 
use the condensed notation for the Green's function of 
the operator 



□ - g^^^x^o 

acting on any symmetric tensor field as 
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with G"^{x,y) - the two-point kernel of this Green's 
function. 

The boundary conditions for this Green's function will 
be discussed in much detail in Sects. III-IV. Here we only 
say that the action ([T]) is formulated in the Euclidean sig- 
nature spacetime, and in the flat-space background setup 
it is understood as metric perturbation expansion on this 
background with trivial zero boundary conditions at in- 
finity which uniquely fix the zeroth order Green's func- 
tion !/□. Formal resummation of this expansion series 
allows one to go over from the fiat-space expansion to 
the expansion on maximally symmetric de Sitter (dS) or 
Anti-de Sitter (AdS) background with the Green's func- 
tion uniquely defined either by the condition of regular- 
ity on S*^ (dS case) or boundary conditions for the Eu- 
clidean AdS spacetime corresponding to the definition of 



^ We use the Euclidean signature spacetime and curvature tensor 



conventions, R = gf'ij^ 
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the Hartle-Hawking vacuum. As the question of our pri- 
mary interest here will be the existence of these two vacua 
(flat-space vacuum and (A)dS one), this information will 
be sufficient to specify the Green's function, for which we 
will only require the following symmetric variational law 
(with respect to local metric variations in □ and P) 
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characteristic of the Euclidean signature d'Alembertian 
with zero boundary conditions. The relevance of this Eu- 
clidean space setup to real physical setting in Lorentzian 
spacetime and, in particular, to causal effective equations 
with retarded nonlocalities will be discussed in Sect. IV 
below. 

The action ([1} has one dimensional parameter M and 
depends on six dimensionless parameters a, a, b, c, d and 
e, the first one a determining the overall magnitude of 
the nonlocal correction to the Einstein term and the rest 
- entering the potential term ^ of the operator □ -I- P. 
For a small value of a and the value of M related to the 
Planck mass Mp, 



\a\ < 1 
Af2 



Ml 



1 



(6) 
(7) 



the theory ([T]) has a GR limit on a flat-space background, 
whereas the rest of the parameters should be restricted 
by the requirement of a stable (A)dS solution existing in 
this theory, which in the dS case can serve as a DE model 
of accelerating cosmic expansion. 

The main claim of our paper is that the theory ([T]) 
has a stable (ghost-free) de Sitter (A > 0) or Anti-de 
Sitter (A < 0) solution of its variational equations of 
motion with an arbitrary A, provided its dimensionless 
parameters satisfy the following restrictions 
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where the new set of capitalized quantities A, B and C 
read as 
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These quantities arise in the coefficients of two tensor 
projectors on traceful and traceless subspaces, which re- 
main in the potential term ([2]), 
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after its calculation on the (A)dS background with the 
Riemann and Ricci tensors of the form 



^ A , . 
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The condition ([5]) guarantees the existence of the 
(A)dS solution, while Eq.(j9l) is responsible for its stabil- 
ity. This fact, which can be regarded as a major technical 
achievement of this work, is based on the calculation of 
the quadratic part of the action on its (A)dS background. 
Bearing in mind generality and complexity of the orig- 
inal multi-parameter action, the result turns out to be 
remarkably simple. When the metric perturbations hfj,i, 
on this background are subject to the background covari- 
ant DeWitt gauge condition. 
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g^^hfj^i, and the indices of h^^ are raised by the 



background metric) this quadratic part depends only on 
the traceless part of hai, and reads 
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Here the condition © signifies the absence of the nonlo- 
cal term and guarantees that the characteristic equation 
for all propagating modes will have only one "massless" 
root for □, □ = 2A/3 (on-shell condition of masslessness 
for (A)dS background differs from its fiat-space analogue, 
□ = [12] )■ Therefore longitudinal and trace modes 
which formally have a ghost nature become unphysical 
and can be eliminated by residual gauge transformations 
preserving the gauge ([TB)l. This well-known mechanism 
leaves us with two transverse-traceless physical modes 
propagating on the (A)dS background - absolutely the 
same situation as in the Einstein theory with two gravi- 
ton modes on a fiat or (A)dS spacetime. 

What is critically different from the GR phase of the 
theory - the effective Planck mass Moff in (fTT]) . which 
determines the cutoff scale of perturbation theory in the 
(A)dS phase and the strength of the gravitational inter- 
action of matter sources. It is given by 
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In the presence of matter sources with the stress tensor 
T^^ the theory on the (A)dS background has linearized 
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gravitational potentials which equal modulo the gauge 
transformation 



S^Geff / □-2A A 

= m , 2A + n , OA ^ ^ ) ■ (20) 



-□- 



□ + 2A SD 



Here all nonlocal operations should be understood with 
retarded boundary conditions (see discussion below in 
Sects. IV and VI) and GcS = l/STrAf^^ is the ef- 
fective gravitational constant vs the Newton constant 
Gjv = l/87rM|,. 

For a wide range of free parameters of our model 
Gcff can be much larger than Gat because in view of 
([5]) a natural range of the parameter B is B ^ a, and 
Geff '■^ Gn/c( ^ Gn- This property can be interpreted 
as a simulation of DM mechanism, because it implies 
strengthening of the gravitational attraction in the (A)dS 
phase of the theory and its possible effect on rotation 
curves at relevant distance scales. 



III. SCALE DEPENDENT COUPLING AND 
FLAT-SPACE BACKGROUND SETUP 

The choice of the action ^ might look contrived, but 
we will show now that it naturally arises within the con- 
cept of the effective scale-dependent gravitational con- 
stant. At a qualitative level this concept was introduced 
in as an implementation of the idea that the effec- 
tive cosmological constant in modern cosmology is very 
small not because the vacuum energy of quantum fields 
is so small, but rather because it gravitates too little. 
This degravitation is possible if the effective gravitational 
coupling constant depends on the momentum scale and 
becomes small for fields nearly homogeneous at the hori- 
zon scale. Naive replacement of the Newton constant by 
a nonlocal operator suggested in violates diffeomor- 
phism invariance, but this procedure can be done covari- 
antly due to the following observation 
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The Einstein action in the vicinity of a flat-space back- 
ground can be rewritten in the form 

where l/D is the Green's function of the covariant 
d'Alembertian acting on a symmetric tensor. This ex- 
pression is nothing but a generally covariant version of 
the quadratic part of the Einstein action in metric pertur- 
bations h^i, on a flat-space background. When rewritten 
in terms of the Ricchi tensor i?^^ ~ Wh + 0[h^] this ex- 
pression becomes nonlocal but preserves diffeomorphism 
invariance to all orders of its curvature expansion. In its 
turn the quadratic nature of the Einstein action in the 
vicinity of a flat-space background follows from the sub- 
traction of the linear in h/^^i, term by the surface Gibbons- 



Hawking integral over asymptotically-flat infinity 

+ — £ / da^^i^d^h^^-d^h). (22) 

From the viewpoint of the metric in the interior of space- 
time this surface term is a topological invariant de- 
pending only on the asymptotic behavior g"^^ = S^^i, + 
hfiuix) I \x\~^oo- It can be converted into the form of the 
volume integral and covariantly expanded in powers of 
the curvature. This expansion starts with^ 



= J dx<?i/2|^-i?^''iG^. + 0[i?;jj|,(23) 



so that the Ricci scalar term gets canceled in ([^^ and 
we come to ([2T|) . 

With this new representation of the Einstein action, 
the idea of a nonlocal scale dependent Planck mass 
can be realized as the replacement of Mp in pT|) by a 
nonlocal operator - a function M^(n) of □, 



S-^MBg,.. (24) 



which would realize this idea at least within the lowest 
order of the covariant curvature expansion. This mod- 
ification put forward in [l5l | did not, however, find 
interesting applications because it has left unanswered a 
critical question - is this construction free of ghost in- 
stabilities for any nontrivial choice of M^(n)? Here we 
try to fill up this omission and put some constraints on 

A'P{a). 

To begin with, if we adopt this strategy, then the search 
for M^(n) should be encompassed by the correspondence 
principle according to which nonlocal terms of the action 
should form a correction to the Einstein Lagrangian aris- 
ing via the replacement R ^ R + R^"'F{D)Gi,^. The 
nonlocal form factor of this correction F{0) should be 
small in the GR domain, but it must considerably mod- 
ify dynamics at the DE scale. Motivated by customary 
spectral representations for nonlocal quantities like 



F{a) = J dm' 



a{m?) 
m? — □ 



(25) 



we might try the following ansatz, F{U\) = a/{m? — □), 
corresponding to the situation when the spectral density 



^ Covariant way to check this relation is to calculate the metric 
variation of this integral and show that its integrand is the total 
divergence which yields the surface term of the above type linear 
in 5g^^{x) = h^,/{x) 
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a{m'^) is sharply peaked around some m? (cf. a simi- 
lar discussion in As we will see, for ^ this 
immediately leads to a serious difficulty. Schematically 
the inverse propagator of the theory - the kernel of the 
quadratic part of the action in metric perturbations hf^i, 
- becomes 



□ 2 



□ ' 



(26) 



where the squared d'Alembertian follows from four 
derivatives contained in the term bilinear in curvatures. 
Then its physical modes are given by the two roots of this 
expression - the solutions of the corresponding quadratic 



equation □ 



In addition to the massless graviton 



with ~ massive modes with — 0{m^) appear 
and contribute a set of ghosts which cannot be eradicated 
by gauge transformations (for the latter have to be ex- 
pended on cancelation of ghosts in the massless sector 
- longitudinal and trace components of hf^i, subject to 

□V = o-)- 

Therefore, only the case of = remains, and as 
a first step to the nonlocal gravity we will consider the 
action 



-R 



(27) 



(for brevity we omit the surface integral that should ac- 
company the Einstein Ricci scalar term). On the flat- 
space background this theory differs little from GR pro- 
vided the dimensionless parameter a is small, \a\ ^ 1. 
Upper bound on jal should follow from post-Newtonian 
corrections in this model. The additional effect of a is 
a small renormalization of the effective Planck mass. In 
the linearized theory we have an obvious relation 



S=- 



dxg^/^R + aO[h 



(28) 



which allows one to relate the constant M to Mp by 
Eq.dll). 



IV. TREATMENT OF NONLOCALITY: 
SCHWINGER-KELDYSH TECHNIQUE VS 
EUCLIDEAN FIELD THEORY 

At this point we have to address the treatment of non- 
locality in (P71) and (P). In principle, handling the the- 
ories having a nonlocal action at the fundamental level 
is a sophisticated and very often an open issue, because 
their nonlocal equations of motion demand special care 
in setting their boundary value problem. Contrary to lo- 
cal field theories subject to a clear Cauchy problem setup 
and local canonical commutation relations, nonlocal the- 
ories can have very ambiguous rules which are critical for 
physical predictions. In particular, the action (j27p above 
requires specification of boundary conditions for the non- 
local Green's function l/D which will necessarily violate 



causality in variational equations of motion for this ac- 
tion. Indeed, the action (|27l) effectively symmetrizes the 
kernel of the Green's function G{x,y) of l/D, so that 
nonlocal terms in equations of motion 



SS 
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{R{y) denoting a collection of curvatures) never have re- 
tarded nature even when G{x, y) is the retarded propa- 
gator or satisfies any other type of boundary conditions 
|10| . Therefore, these equations break causality because 
the behavior of the field at the point x is not independent 
of the field values at the points y belonging to the future 
light cone of x, y° > 

To avoid these ambiguities and potential inconsisten- 
cies we will once and for all assume that our nonlocal 
action is not fundamental. Rather it is the quantum ef- 
fective action - the generating functional of one-particle 
irreducible diagrams - whose functional argument is the 
mean quantum field. This functional is necessarily non- 
local, and its nonlocality originates from quantum effects 
(by various mechanisms widely discussed in literature in- 
cluding [23, UHl)- In this case boundary conditions for 
nonlocal operations are uniquely fixed by the choice of 
the initial (and/or final) quantum state, and manifest 
breakdown of causality in variational equations for this 
action is harmless under a proper treatment of their non- 
local terms. 

To begin with, this causality breakdown does not im- 
mediately signify inconsistency in the calculation of scat- 
tering amplitudes or ( in | out ) matrix elements, because 
these amplitudes are determined by Feynman diagram- 
matic technique and do not have manifest retardation 
properties because they are not directly physically ob- 
servable. Physically observable quantities like probabil- 
ities are bilinear combinations of scattering amplitudes 
and can always be represented as expectation values 
( in I O I in ) of certain quantum operators O in the ini- 
tial quantum state | in). For example, the probability of 
transition from this state to some final state I fin ) 



Pi, 



>fin = { m 



fin )( fin I 



= ( in I p. 



fin I m; 



(30) 



is an expectation value of the projector Pgn = | fin ) ( fin | 
onto this final state. In contrast to in-out matrix ele- 
ments these expectation values are subject to Schwinger- 
Keldysh diagrammatic technique (l9l . [20| which guaran- 
tees causality of (in | ©(a;) | in). This property can be 
formulated as a retarded response of this average to the 
variation of the classical external source J{y) coupled to 
the quantum fields in terms of which the operator 0{x) 
is built, 



5{in\d{x)\ 



0, 



(31) 



This property is again not manifest and turns out to 
be the consequence of locality and unitarity of the origi- 
nal fundamental field theory (achieved via a complex set 
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of cancellations between nonlocal terms with chronologi- 
cal and anti-chronological boundary conditions). How- 
ever, there exists a class of problems for which a re- 
tarded nature of effective equations of motion explicitly 
follows from their quantum effective action calculated in 
Euclidean spacetime [2l|. This is a statement based on 
Schwinger-Keldysh technique that for an appropri- 
ately defined initial quantum state |in) the effective equa- 
tions for the mean field 



^mlg^^lm; 



(32) 



originate from the Euclidean quantum effective action 
S = ^EuciidcaniS/ii/] by the following procedure [2lj^. Cal- 
culate nonlocal SEucmc&nig fii^] and its variational deriva- 
tive. In the Euclidean signature spacetime nonlocal 
quantities, relevant Green's functions and their varia- 
tions are generally uniquely determined by their trivial 
(zero) boundary conditions at infinity, so that this vari- 
ational derivative is unambiguous in Euclidean theory. 
Then make a transition to the Lorentzian signature and 
impose the retarded boundary conditions on the resulting 
nonlocal operators. 



(55, 



Euclidc 



retarded 



0. 



(33) 



These equations are causal {g^v{x) depending only on 
the field behavior in the past of the point x) and satisfy 
all local gauge and diffeomorphism symmetries encoded 
in the original S'Euciidcan[5pi']- 

A similar treatment of a nonlocal action in [lo| was 
very reservedly called the "integration by parts trick" 
needing justification from the Schwinger-Keldysh tech- 
nique. However, this technique only provides the causal- 
ity of effective equations, but does not guarantee the 
Euclidean-Lorentzian relation ((33|) . The latter is based, 
among other things, on the choice of the |in)-state. 

We will assume that our model falls into the range 
of validity of this procedure, which implies a particular 
vacuum state |in) and the one- loop approximation (in 
which it was proven to the first order of perturbation 
theory in [26| and to all orders in [2l| ) . The extension of 
this range is likely to include multi-loop orders and the 
I in) -state on the (A)dS background considered below, for 
which this state apparently coincides with the Euclidean 
Bunch-Davies vacuum. 



crating DE. Though it will fail due to inconsistent treat- 
ment of boundary conditions, it is instructive to pass this 
exercise to see the importance of their careful treatment. 

Like in papers on /(i?/n)-gravity (see [l^] and refer- 
ences therein) stemming from Q, but in contrast to Q 
disregarding consistent treatment of boundary conditions 
for nonlocal operations, one can localize the nonlocal part 
of (|27p with the aid of an auxiliary tensor field (^^^ . Then, 
the theory is equivalently described by the action 

-a(^'''^-i5^>)n^^.} (34) 

generating for ip^^'^ the equation of motion □(^'^'^ — —G'^'^. 
Since ([27]) is understood as the Euclidean action with zero 
boundary conditions for !/□ at infinity, the auxiliary 
tensor field should satisfy the same Dirichlet boundary 
conditions (^^'^|oo = 0, and this is critically important 
for stability of the theory. Indeed, the field ip^'^ formally 
contains ghosts, but they do not indicate physical in- 
stability because they never exist as a free fields in the 
external lines of Feynman graphs. In the Lorentzian con- 
text of ([55]) this means that yj^" is given by a retarded 
solution, ip^" = —{l/Dj^ctG^" , and does not include free 
waves coming from asymptotic infinity. 

Artificial nature of these ghosts is analogous to the 
case of the simplest ghost-free action that can be formally 
rendered nonlocal 



dx (pOip 



dx (□<^)i(n^) 



and further localized in terms of the auxiliary field "0 with 
the action 

S'[ ^, ?A ] = J dx (21/- Dtp + 11! □V') . 

This action is equivalent to the original one when is 
integrated out with the boundary conditions (ip+tp) \ oo — 
0. After diagonalization this action features the ghost 
field g = ip + ip, 



S[ip,^] — J dx (g Dg — ip Dip) . 



(35) 



A. Asymptotically-flat space vs cosmological setup 

This subsection demonstrates trial application of the 
model (j27p in cosmological setup with the purpose of gen- 



^ We formulate this statement directly for the case of gravity 
theory with the expectation value of the metric field operator 
Qiiui^)-, though it is valid in a much wider context of a generic 
local field theory [S^]. 



This ghost is however harmless because under the bound- 
ary conditions of the above type it identically vanishes in 
view of its equation of motion Dg — 0. In the presence 
of interaction, a nonvanishing g exists in the intermedi- 
ate states, but never arises in the asymptotic states, or 
external lines of Feynman graphs. 

Main lesson to be drawn from the above example is 
that the actual particle content of the theory should be 
determined in terms of the original set of fields, whereas 
nonlocal reparameterizations can lead to artificial ghost 
modes which are actually eliminated by correct boundary 
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conditions.^ In our case this is the original formulation 
(P7)) in terms of the metric field g^^. It indeed turns out 
to be ghost-free on the flat-space background, because 
the quadratic part of the action coincides with the Ein- 
stein one.^ 

In the local representation (p4|) our model could be di- 
rectly applied to the FRW cosmology for the purpose of 
finding the accelerated stage of cosmic expansion. It is 
easy to find a (quasi) de Sitter point of the cosmologi- 
cal evolution. Indeed, with the natural Lorentz-invariant 
ansatz for the auxiliary field (p'^'^ ~ g^"^ , which is sup- 
posed to be valid close to a certain moment to correspond- 
ing to the present epoch, the cosmological evolution for 
the action can be compatible with the current DE 
data. By an appropriate choice of initial conditions the 
Hubble factor H = a /a, the field $ and the parameter 
of the effective equation of state w = — 1 — 2H/3H^ can 
satisiy at to the following relations 



a 



Wo 



'1, wo = -l&Ho 



2cr - 1 

30-2 + 2 



(36) 



Here the quantity a = {2a/ 3)^^'^ {2 + a<^>o - 3a) "^/^ jg 
determined by the value of the field $0 = ^{to)- If it is 
chosen to satisfy a = 0(1) > 1/2 we have Wo — 0(1) x 
Hq < which yields a pace of change in the effective 
equation of state compatible with the horizon scale. 

These preliminary estimates could have served as a 
starting point for a quantitative comparison with the DE 
scenario. However, a formal application of (P7|) to the 
FRW setup disregards nontrivial boundary conditions in 
cosmology. To see this, note that fine tuning initial con- 
ditions for $ to the DE data would generally contra- 
dict zero boundary conditions for the auxiliary tensor 
field (fi^'^ cx ^g^'^, not to mention that the cosmolog- 
ical FRW setup does not in principle match with the 
asymptotically-fiat framework of the action (P7)) . There- 
fore we have to extrapolate the definition (l?7l) to non- 
trivial backgrounds including, first of all, the de Sitter 
spacetime and change our technique - instead of localiza- 
tion method with an auxiliary tensor field work directly 
in the original metric representation. This will help us 



This, in particular, means that the ghost avoidance criteria de- 
rived in I27II are not precise, because they are in fact applied to 
local scalar-tensor theories rather than to nonlocal ones. We have 
to reiterate here that a naive analysis of kinetic terms of auxiliary 
fields, which are usually used to localize a nonlocal action of the 
theory, cannot excusively serve as a criterion of the elimination 
of ghosts — boundary conditions are equally important for that. 
A similar mechanism of eliminating ghosts by boundary con- 
ditions was recently used in [28l ] . However, in contrast to our 
model the ghost modes of the conformal gravity in [28l | are higher- 
derivative ones and are essentially nonlinear. Therefore, the non- 
ghost nature of the theory requires further verification even after 
integrating these ghosts out . 



to look for the (A)dS solution in the covariant language 
circumventing explicit FRW metric ansatz. 

This approach to the action (|27p suffers from a seri- 
ous difficulty. Ricci curvature for the (A)dS background 
(|15p is covariantly constant, and the nonlocal part of 
([27| turns out to be infrared divergent, {l/0)gf^^ — 00 
- the property that can hardly be cured by some choice 
of boundary conditions for l/D. This means that the 
action (|27|) should be modified to regulate this type of 
divergences. 



V. STABLE (A)DS BACKGROUND 

We will regulate the action (|27|) by adding to the co- 
variant d'Alembertian the matrix-valued potential term 
built of a generic combination of tensor structures linear 
in the curvature. This brings us to the six-parameter 
family of nonlocal action functionals introduced 
in Sect. II. Of course, such a modification of the origi- 
nal action (|27l) leaves its linear approximation on a flat 
background intact, because it deals with 0[/i^^]-terms, 
whereas its dimensionless parameters will be restricted 
by the requirement of a stable dS or AdS solution in the 
model. 

The action of the matrix valued potential P on g^^, is 
given by the relation 

Pg^, = P^^''^ g^0^AR^, + BRg^,, (37) 

where the coefficients A and B are defined by Eqs. dTU]) - 
dTTI) . so that the Green's function 1/{D + P) acting on 
Ricci and Einstein tensors in ([T]) is well defined even for 
the maximally symmetric (A)dS background with the co- 
variantly constant curvatures p^ -([T5 |) . In this case the 
above relation simplifies to the equation 



Pgf,^ = {A + 4B)Ag^^ 
which has two obvious corollaries 



1 



D + P 
1 

D + P 



Rn 



AP- 



1 



(A)dS 



(A)dS 



9^a 



1 



A + AB 



(38) 



A + AB 



(40) 



These corollaries follow from the fact that on a maxi- 
mally symmetric background □ commutes with P and 
annihilates i?^^ = Ag^^. 

Let us prove now that under the relation ^ between 
the parameters of the potential ([2]) the equation of motion 
for ([l} has the (A)dS solution with an arbitrary value of 
the cosmological constant A. The simplest method to 
see this is to apply to the action the conformal metric 
variation (more cumbersome straightforward calculation 
will be presented in Appendix). 

For this introduce the local conformal variation with 
the parameter 5a — Sa{x), 

d'^xSagajS^ ■ (41) 

Ogal3 
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Under its action various quantities in ([l} transform ac- 
cording to their conformal weights, 

S^Rf"" = -2(5ct ff'" + 0(V), 5aP = -6aP + 0(V), 
d^n = 0(V), (42) 

modulo the derivatives 0(V) acting on ^17(0;) and these 
quantities themselves. Then the conformal variation of 
([T]) on the (A)dS background reads 

S,S=^ J d^xg^^^ (^-R + aR''^-^^Gc,p\5a 



= -2M^A 1 



A + AB 



^ J d^xg^^^Scr, (43) 



so that the corresponding variational derivative equals 



SS_ 



= -2Af^A 1 



(A)dS 



A + 4B 



,1/2 



(44) 



Since all tensor quantities on this background alge- 
braically express via the metric variational derivative 
of the action reduces to this conformal variation. Then 
the resulting equation of motion 



(55 



(A)dS 



4^ (5a 







(45) 



(A)dS 



holds with an arbitrary value of the effective cosmological 
constant A when the parameters in Eqs.([T])-([2|) satisfy the 
relation ([S]). 

Note that the existence of the (A)dS solution with an 
arbitrary A is neither the result of the local Weyl invari- 
ance of the theory, nor even its global scale invariance. 
Rather this is the corollary of the relation ([5|) which, in 
virtue of Ea.(l40)). guarantees the vanishing on-shell value 
of the action. 



^ I (A)dS ^ ^ 



(46) 



Thus, this solution is another vacuum - a direct analogue 
of the flat-space one. 

Another remarkable corollary of Eq.Q is that the sta- 
bility of the (A)dS solution against ghost and tachyon ex- 
citations is guaranteed by only one additional restriction 
© on the parameter C defined by p^ . In principle, the 
hope to eliminate ghosts and tachyons from the quadratic 
part of the action S(2) on the (A)dS background is based 
on the following observation. 

In the DeWitt gauge on metric perturbations 5(2) 
contains only two structures /I'^^x /i^^ and /ix /i, or equiv- 
alently h^'^ x h^^y and hx h, where /i^^ is a traceless part 
of hfj^^. Since the potential term P on the maximally 
symmetric background commutes with □ and equals the 
linear combination (jl3p of projectors on traceless and 



trace subspaces, the nonlocal parts of these two struc- 
tures have the form 



□ -CA 



D~aA 



(47) 



where the Green's function in the trace sector follows 
from the equation 



[D + P) g^,h ^ g^,{a - aA) h 



(48) 



which is also based on dHj) and (p8l) . As in the flat-space 
background discussion above, the ghosts necessarily ap- 
pear if these nonlocalities are nonvanishing in 5(2), be- 
cause the dispersion equation for □ becomes quadratic 
and generates a doubled set of physical modes with 
□ = TOj_ (with CA or a A playing the role of a nonvan- 
ishing m^). Therefore it seems a priori possible to cancel 
these two nonlocal structures and provide the right signs 
of the remaining local terms by the appropriate choice of 
the five free parameters in 

Curious fact is that in the DeWitt gauge 5(2) has 
a very simple form 



5(2) = ^ j d^xg^/^i^--^h^^'^Uh^, + ^^hUh 
-i (c~^\Ahl + ^(c-^\Ah^ 



(49) 



U + P 
4 U-aA 



where the effective Planck mass M^ff is given by Eq. ([T9| . 

A rather lengthy technical derivation of this result is 
presented in Appendix. It is strongly based on the re- 
lation ([8]), Eqs. (p9l) -(|40 |) and, what is less pronounced 
but equally important, on a symmetric variational rela- 
tion for the Green's function ([5]) and the possibility of 
integrating by parts without extra surface terms. These 
properties are guaranteed by the Euclidean spacetime sig- 
nature in the action ([T|) and regularity of Green's func- 
tions on a closed compact S'^ for the de Sitter case or 
their boundary conditions at the AdS boundary for the 
Anti-de Sitter case. 

Again, it is worth mentioning here that this rule of free 
integration by parts, that was used throughout the series 
of papers fol-lll | and interpreted there as a trick, makes 
sense only within the relation between the Schwinger- 
Keldysh technique and Euclidean QFT and applies only 
in the Euclidean spacetime. Below we will see that, once 
the effective equations in the physical Lorentzian space- 
time have been obtained by Eq. psp , indiscriminate omis- 
sion of surface terms under integration by parts becomes 
illegitimate and leads to wrong properties of free propa- 
gating modes of h^i, . 

Remarkably, the expression (|49p depends only on the 
traceless part (ITS)) of the metric perturbation, because 
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the combination of tensor contractions hi^" x /i^^ — jhxh 
equals h^'^ x hfj^^. This property, apart from expUcit 
calculation done in the Appendix, can be derived from 
the scaling transformation of the action ([T]) 



5[(l + £)g^,.] = (l+e)5[g^,] 



(50) 



under a global dilatation of the metric, V^e = 0. This 
implies that the quadratic part of the action (in fact on 
an arbitrary background) identically vanishes for global 
conformal perturbations h^i, = e gf_iu obviously satisfying 
the DeWitt gauge ^ 



(2)1 



0, 



const. 



(51) 



This observation serves as an independent check of ex- 
plicit calculations of 5'(2) in Appendix and uniquely de- 
termines the minus one quarter ratio of the coefficients 
oi h X h and /i^" x /i^i, structures (except their local 
kinetic terms hOh and h'^'^Dhf^^ which identically vanish 
for h^i, — £.9^1/ with a constant e). 

Thus in the DeWitt gauge the quadratic part of the 
action takes the form (IT71) quadratic in the traceless 
part of metric perturbations. Moreover, in view of the 
projector properties of the potential term (I13p it simpli- 
fies even further 



Ml 



(2) 



d X g 



□ -CA 



(52) 



With an arbitrary value of C the characteristic equation 
for physical modes still yields only one root □ = |A, but 
this is the second order root which for C 7^ | corresponds 
to double poles in the propagator and signifies presence 
of higher-derivative ghosts. 

Therefore, the requirement of absence of ghosts im- 
poses only one extra equation ^ for C, C = 2/3, and 
the positivity requirement for M'^^. Bearing in mind that 
I a I <C 1, this selects two admissible intervals for the pa- 
rameter B in the case of a positive a, 



a 



B<--, B>0, 



(53) 



and even more interesting compact range of this param- 
eter for a negative a. 



0<B < , a < 0. 

2 



(54) 



VI. FREE PROPAGATING MODES AND 
GRAVITATIONAL POTENTIALS IN THE (A)DS 
PHASE 

In order to obtain correct linearized equations of mo- 
tion we need a quadratic part of the action without 
gauge-fixing. It is invariant under the linearized gauge 
transformations h^^i, — >■ /i^^ -I- A^/i^^ with the vector dif- 
feomorphism parameters = g^vf^ , 



(55) 



where covariant derivatives are determined with respect 
to the background metric. The invariant S'(2) can be 
obtained from by representing h^'^ in the DeWitt 
gauge as the projection to this gauge of the non-gauged 
field. 



2V 



x°=o 



□ + A 



Xo 



(56) 



Here it is worth reminding the definition of the DeWitt 
gauge condition functions (jl6l) for which the vector field 
operator (□ -I- A) plays the role of the Faddeev-Popov 
operator 



1 



A/x'^ = (n + A) r. 



(57) 
(58) 



The result of substituting in Eq.(|49l) for C 



reads 



S'(2) - 



:,h[-U---K] h--xl 



4^(1) 



1 



□ + 2A 



^(1) 



(59) 



where is the linearized Ricci scalar on the (A)dS 
background which has a form of the combination of two 
terms linear respectively in V^x'^ and h, 

- V^X^ - ^ (° + 2A) h. (60) 

Interestingly, the first two lines of S'(2) above coincide 
with the quadratic part on the (A)dS background of the 
Einstein-Hilbert action with the A term [22|- Therefore, 
this part of 5(2) is invariant under (j55p . while the invari- 
ance of the last term of ((59)) directly follows from the 
invariance of R(i)- The form of ([5^ shows that despite 
a complicated tensor structure of nonlocal term in ^ 
near the (A)dS background it reduces to the Ricci scalar 
sector similar to the R-^R distortion of Einstein theory 
considered in . 

Now, according to the relation ([55)) between Euclidean 
QFT and causal effective equations for mean field we cal- 
culate the variational derivative of S'(2) with respect to 
/i^^, go over to the Lorentzian metric signature and es- 
tablish retardation of all nonlocal operations. This gives 
the following equation for free propagating modes of the 



( expectation value ) metric field h 



•■fj.i' I 



-1/2 '^'^(2) 



retarded 



h+++ -+++ 



1 



-I- 2 V(^<Z',) - g^.V„<?" = 0, (61) 



10 



where is the nonlocal function 

1 



1 v^- 

2 □ 



2A 



(62) 



whose nonlocality is given by the retarded Green's func- 
tion . 

The trace of this effective equation gives 



Mis 



9'^ 



□ 



1 



^(1) - 



2A 
2A 



R(i 



(1) 



□ + 2A 



- 0. (63) 



The last equation yields not only the homogeneous dif- 
ferential equation for but also its zero initial condi- 
tions at past infinity because of the retarded nature of the 
Green's function. Acting on this equation by the opera- 
tor □ + 2A we immediately get the initial value problem, 



= 0, 



(64) 

%)Lo^-oo = 0' ^oi?(i)Lo^_=0, (65) 



with the identically vanishing solution. Therefore the 
linearized Ricci scalar of the free propagating wave is 
vanishing throughout the entire (A)dS spacetime 



Rii){x)=0. 



(66) 



pairs of initial data for accomplishes the counting of 
the physical degrees of freedom among spatial compo- 
nents of hf^^, 6—1 — 3 = 2, while the four lapse and 
shift functions /iq^, as usual, express via the constraint 
equations of motion SS(2)/Sho^ — 0.^ 

It is important that a priori in the DeWitt gauge the 
equation of motion (|6ip is not local even though the func- 
tion (p^ seems getting localized, cf. Eg . ([SO)) . This is 
because the seemingly correct equation 



1 



□ + 2A 



( 1^ + 2A) (^(.t) = ip{x) 



(70) 



does not hold for the field (p{x) with a noncompact sup- 
port extending to past infinity. This is exactly the case 
of free propagating modes, when integration by parts 
of □ yields extra surface terms - the situation differ- 
ent from the Euclidean QFT side of the relation (|33|) . 
as it was mentioned in the previous section. Disregard- 
ing these terms would lead to inconsistent equations of 
motion breaking their diffeomorphism invariance.^ 

Opposite situation occurs when we consider retarded 
gravitational potentials /i^i, generated by matter sources 
with a stress tensor T^^ having a compact support. Then 
/i^ii/ also has a compact support in the past, and Eg. ((70)) 
applies to the calculation of the vector function (|62|) . In 
the DeWitt gauge it becomes local, <?''|xa=o = -jV^/i, 
and the equation for the retarded potential also takes 
the local form. 



As a result the nonlocal function ([62)) coincides with the 
local DeWitt gauge condition function, = x^, and 
Eq. (|6T)) becomes absolutely identical to the linearized 
Einstein equations on the (A)dS background. 

Now it is high time to impose the DeWitt gauge (fT6)) 
in which Ea. ljMl) reduces to 



(67) 



Similarly to the Feynman gauge in electrodynamics, in 
this relativistic gauge all components of h^^ are propa- 
gating, but their gauge ambiguity is not completely fixed 
and admits residual gauge transformations (|55p with the 
parameter satisfying the equation 



(□-f A)/^-0. 



(68) 



By the usual procedure these transformations can be 
used to select two polarizations - free physical modes 
f^fil^'' = ^M'' + '^tJ.fiy + ^vfp.- In particular, they can 
nullify boundary conditions for h on any initial Cauchy 
surface E, so that this trace identically vanishes in view 
of the homogeneous equation ()67l) and makes the physical 
modes transverse and traceless as in the Einstein theory 
with a A-term, 



Vkfj"" = 0, /jp^y^ = 0. 



(69) 



Indeed, under these transformations hP"^^^ /i + 2Vp/'^, 
and both V^/'^|s and 9oV^/^|s can be chosen to pro- 
vide zero initial data for h^^^^ on E. The remaining three 



-□ + -A 1 /i 



1 A ^ _ 2 



■T^.. (71) 



cff 



with matter stress tensor coupled to hf^i, via the effective 
Planck mass (fT9)). 

The retarded solution of this equation for the trace 
part of hfj,y immediately reads 



T. 



(72) 



After careful commutation of covariant derivatives and 
the Green's function of the tensor operator (-□-!- 



^A)(5 

3 -"^/"/ijy I 



5 



2A 

3 



v„v 



2A 



□ 



-9iJ^u- 



-□ + |A) (□ + 2A)' 



(73) 



" Or equivalently, when they are treated as propagating modes 
subject to second order in time differential equations, their initial 
data express via x'^l E = a^nd t^oX^I E = 0. 

Indiscriminate use of Eq. jTOj, which leads to the local equation 
Il61|l in the DeWitt gauge, would imply the wrong trace equation 
□h = different from H67| l. This would contradict the corollary 
of Eq.imil, (□ + 2A)V^/f = 0, because h and feP^ys differing 
by 2V^/'^ would not satisfy one and the same equation as they 
should. 
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one finally finds the gravitational potential of matter 
sources. In the DeWitt gauge it takes the form 



167rG, 



cff 



3 



□ -2A A 



(□ + 2A) □ 



□ + 2A SD 



T. 



(74) 



Here of course retardation prescription applies to all 
Green's functions, the last term represents a pure gauge 
transformation and Gcs = 1 /^TrM^g is the effective grav- 
itational constant vs the Newton one Gn = l/SnMp, 



GcB — 



a{l — a) 
8B{2B + a 



Gn- 



(75) 



This proves the gauge invariant part of this potential 
advocated in Sect. II. 

The expression ([71)) differs from the GR analog by the 
tensor structure - the nonlocal combination in the first 
line of (IHl) vs the GR structure T^,, - ^gtivT. From 
the viewpoint of expansion in \a\ <^ 1 this implies for 
non-relativistic sources certain 0(1) corrections. What is 
much more interesting, it yields an unexpected bonus in 
the form of a possible dark matter simulation - 0(1/ jal) 
amplification of the gravitational attraction due to the 
replacement of the Newton gravitational constant Gn 
by Goff ~ Gjv/|q;| with |a| ^ 1. This necessarily hap- 
pens in the case ([M]) of a negative a, because the factor 
a/8B{2B + a)> l/\a\ and 



G. 



cff 



> 



1 



GAr>G 



a 



N- 



(76) 



For a positive a in the domain ([55)1 the theory is also 
strongly coupled with Gog > Gn for the following two 
intervals of B, 



<B < 



a 



0<B < 



la - 



(77) 



and might correspond to the DM phenomenon. In the 
rest of this domain it is on the contrary weakly coupled 
with Geff < Gn- At the crossover between these two re- 
gions of strong and weak coupling, \B\ ~ -\/a/4, both 
Newton and effective gravitational coupling constants 
can be of the same order of magnitude, Gn/GcR = 0(1) 
even for a very small a ^ 1, which together with ([7S)) 
leaves a large window for a possible strength of DM at- 
traction relative to the GR behavior. 



VII. CONCLUSIONS 

Thus we have a class of generally covariant nonlocal 
gravity models which have a GR limit and also a stable 
de Sitter or Anti-de Sitter background with an arbitrary 
value of its cosmological constant. This background car- 
ries as free propagating modes two massless gravitons 



identical to those of Einstein theory with a cosmologi- 
cal term. These gravitons are coupled to matter with 
a variable effective gravitational coupling which can be 
stronger than the Newton coupling in general relativity 
theory. The interpretation of these models looks as fol- 
lows. 

The theory with the action ([1]) has two phases. For 
short distances corresponding to the range of wavelengths 
with 



VV - □ > i? 



(78) 



this is a GR phase on the zero curvature background with 
small 0(a) x R/U\ corrections of higher orders in space- 
time curvature (collectively denoted by R). This regime 
might apply to galactic. Solar system and other small 
scale phenomena and is likely to pass all general relativis- 
tic tests for a sufhciently small a. Disturbing property 
of perturbation theory in this phase could be the pres- 
ence of poles ~ (l/n)i? in the vertices of ([TJ which could 
contribute too strong to graviton scattering. But these 
contributions vanish on shell R^^ — provided, perhaps, 
that the Riemann term is forbidden in the potential © of 
the tensor differential operator which specifies a nonlocal 
part of (P), a = 0. 

Another phase of the theory corresponds to the in- 
frared wavelengths range 



VV < i? 



(79) 



in which a stable (A)dS background exists and the mod- 
ified gravitational potential of matter sources is given by 
Ea. (l74)) . This equation is valid for the perturbation range 
I (5i?^ I I VV/i(; I < A and I I < 1, which is equiva- 
lent in virtue of Eq. (j7ip to very small matter densities, 



(80) 



characteristic of galaxy, galaxy cluster and horizon scales 
for which DE and DM modification of gravity theory be- 
comes important. Thus, nonlocal gravity interpolates be- 
tween GR theory and its infrared modification. The lat- 
ter is likely to generate a stable ghost-free stage of cosmic 
acceleration and, perhaps, even simulate the DM effect 
on rotation curves in a strong coupling domain (j54p and 

dZZl). 

How realistic is this picture from the viewpoint of the 
concordance model of DE and DM? There are open prob- 
lems which might derail its viability. The most serious 
objection is that absence of ghosts is guaranteed only on 
two backgrounds ~ fiat space and (A)dS ones. Outside 
of these two solutions the theory is most likely unstable, 
which perhaps can be interpreted as a kind of attraction 
mechanism pushing the system to one of its two stable 
phases. Therefore, inclusion of matter sources can be 
consistently done only within perturbation theory in the 
vicinity of these two solutions, that is when the total mat- 
ter stress tensor T^^ is treated as a perturbation. In fact, 
the range of energies (j80p is a condition justifying this 
perturbation theory on the (A)dS background. A similar 
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range on the flat background is of course bounded by the 
Planck scale which is not at all restrictive, because we 
anyway stay in the infrared domain of effective gravity 
theory. 

It goes without saying that the model should undergo 
tests on consistency of post-Newtonian corrections, the 
magnitude of the DM phenomenon should be compared 
with observations, the effect of nonlocal stress tensor 
trace term of ((74|) should be studied, etc. Moreover, the 
mechanism should be found, by which the model picks 
up a concrete scale | Tj^ \ ~ Af^j^A of a crossover from the 
GR regime to cosmic acceleration - necessary element in 
realistic cosmology. Regarding this mechanism we only 
note that the existence of an (A)dS solution with an ar- 
bitrary A is possible due to the fact that the purely grav- 
itational action ([l} transforms homogeneously under the 
global rescaling g^,^ \g^,y, S[\g^,v] = \S[g^,u\-, so 
the crossover mechanism can be based on the breakdown 
of this property by matter fields, cf. 

The last but not the least is the justification of the 
choice of (H]) as an approximation for the effective action 
coming from some fundamental quantum gravity theory. 
Quantum effective action with the scaling property of 
the above type is hard to imagine within semiclassical 
expansion which contains additional to Mj, dimensional 
parameters and has another scaling behavior [U, [25j . 
However, infrared nonlocal expansion in the heat kernel 
theory of grav itating models and their nonlocal effective 
action |23l |29| contain nonlocal structures similar to ([I} . 

Apart from pragmatic applications in cosmology the 
model in question can be interesting from pure field- 
theoretical viewpoint. It is usually considered that mod- 
ifications of Einstein theory are associated with new de- 
grees of freedom of the gravitational field or broad gravi- 
ton resonances [l3l . The above model shows that it is not 
necessarily the case, because both phases of the theory 
have the same set of Einstein massless modes. In fact, 
our model presents the way how the vacuum equation of 
motion 



R^^ - Kg^i, = 



(81) 



can be encapsulated into the action functional in such 
a sophisticated way that feeding the theory with mat- 
ter sources takes place with a nonlocal coupling varying 
from one of its stable phases to another. Usually this is 
very hard to implement if we want to maintain stability 
and unitarity of the theory. Interesting example of such 
a procedure is a local model of conformal gravity [1^ , in 
which higher-derivative ghosts are advocated to be elim- 
inated by special boundary conditions. Our model is the 
example of a nonlocal model explicitly demonstrating ab- 
sence of ghosts on two vacua of the theory and possible 
strong and weak coupling regimes associated with these 
vacua. Moreover, both of these flat and (A)dS vacua have 
zero value of their on-shell action (|46|. which makes the 
situation very attractive because their contributions to 
physical amplitudes become equally important without 
being inflnitely suppressed or enhanced, as it happens in 



other models of Euclidean quantum gravity.^ 

In conclusion we mention that serendipity of ghost-free 
nonlocal gravity models ([T|) might not be exhausted by 
applications in cosmology. In particular, they might ad- 
mit generic (not maximally symmetric) Einstein space 
solutions (|81l) with an arbitrary cosmological constant. 
Therefore, these models can have implications in black 
hole physics and be an alternative to the conformal grav- 
ity model of [11] , whereas with a negative A they become 
a new testing ground for AdS/CFT correspondence per- 
haps promising other exciting consequences. 
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NOTE ADDED IN PROOF 

At the proofreading stage of this work the paper [s^ 
appeared, advocating that the above results can be gen- 
eralized to the case of a generic not maximally symmetric 
Einstein space background, provided that the coefficient 
of the Riemann term a in Eq.(2) equals a = 2 (in contrast 
to the initially expected zero value). This has important 
implications for a wide class of Schwarzschild-de Sitter 
black hole solutions with vanishing horizon entropy {30| . 



Appendix A: Quadratic part of the action 

Symmetric variation rule for the Green's function ([5]) 
allows one to write down the first order variation of the 
action ([T} under the metric variation Sg^^i, = h^^i, 

SS^^ f d^xg^'^\G^''h^, + ^hR^^''-^G^, 

z J L LJ H- _P 

□+F D+P 

-aff"'^^(5^ + 5P)^-.G^,l. (A.l) 

□ + P ^ ' U + P \ 



* The degeneracy in the family of (A)dS vacua with different A 
might be interpreted as instabiUty with respect to classical or 
quantum leaps between them. However, classically A is a con- 
stant of motion, whereas its quantum mechanical treatment does 
not make much sense, because we already work within effective 
theory for mean fields, which resulted from quantization, and A 
is the mean value of the effective cosmological constant. 
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Here and in what follows we use obvious but deserving 
special mentioning notations 



The arrow over □ indicates the direction in which this 
operator is acting. For its metric variation SO it is im- 
portant to indicate whether it acts to the right or to 
the left, because in contrast to the symmetric operator 
g^^^O the operator g^^^SD is not symmetric, and under 
integration by parts (reversing the direction of the arrow) 
generates extra terms cx Sg^^^ = ^g^^^h. 

The expression (jA.ip allows one to derive explicitly the 
variational derivative of the action on the (A)dS back- 
ground and confirm Eas. p4)) - (|45]) derived above by the 
conformal variation method. Using (|39|) and (|40|) and 
reversing the action of the symmetric {O + P)~^ to the 
left in the fourth and fifth terms of (|A.ip we have 



SS 



(A)ds 2 J ^ \ A 



A 



4B 

a 



gp.y 



4B 



[A + 45) 



■g 



A + iB 
(5^ + SP) g^. 



(A.2) 



To calculate the first four terms here we use the Ricci 
tensor variation 



-Oh, 
2 



iv^V./i, (A.3) 



while the fifth term can be obtained as follows. Consider 
the first order metric variation of the identities Og^,^ = 
and ([57| . This gives 



{S~^) g^u = -Oh^j.v 
{5P)g^,^~Ph 



f^i, -r A 5R^u 

B gf_,„ SR, 



BRh. 



so that on the (A)dS background 

{S'^ + SP)g^, ^~0 + P)h 



ASR 



-ABA h^^ + Bg^^ SR. 



(A.4) 
(A.5) 

(A.6) 



Therefore the variation (IA.2I) takes the form 

APK / a 
(A)ds ~ 2~ V ^ A + AB 



SS 



d^xg^'^g^^'^Sg^, 
(A.7) 

and confirms Eqs. (PH) - (H5]) which lead to the criterion ([5]) 
of the (A)dS background with an arbitrary A. 

Now let us go over to the quadratic part of the ac- 
tion on this background. For this we make the metric 
variation of (lA.ip by applying the same rule ([5]) and in- 
tegration by parts 



S{2) = ^s'^s 



(A)dS 



(A., 



The result immediately simplifies if we use the following 
two corollaries 



g- 



1 



U + P 



1 



□ - aA 



(.9'-''$^,), (A.9) 



(A.lO) 



of the equation 

g-P{U + P)J = {U-aK)g^'^, (A.ll) 

which is of course based on the relations and ([5]). 
The above equations hold for generic tensor ^^'^ and 
scalar $ fields and allow us to pull the metric tensor 
through the tensor propagator converting it into the 
scalar one. 

As a result the quadratic part of the action splits into 
the sum 

S(2) = <5'(2) + 5'(2) (A. 12) 

of the purely local term 

S(,)^^ j d'xg''^i^\hSR+\Kh^ 

+5G^ni^, - i hl^R + i hg^JR^-" 

+ {S^R^''')g^,,~g''''{S^G^,) 

+ igf^-{S^^ + S'p)g^,'j (A.13) 

and the term which together with some local pieces con- 
tains all nonlocalities of iS'(2) 

^ [a 1/2 / , «A ^ 1 
S(2) = — r I xg ' <j +—h— —g^ SG^^ 



4 ./ " " " I ' 2 □ - aA ■ 
A , 1 



■g 



((5^ + SP) g^. 



2 □ - aA ■ 



a 



-2aSR'"'- 



5G,. 



□ + P 

-2SR^"'—^ {st\+SP) 



U + P 



gilv 



+2g^'''{S'^ + 5P)^^SG^}^. (A.14) 



1. Disentangling nonlocal structures 

Here we begin with disentangling nonlocal terms from 
(|A.14[) . From now on we will consider metric perturba- 
tions in the DeWitt gauge in which the Ricci tensor 
variation (|A.3p (on the (A)dS-background) simplifies to 



SRiiv 



(^^,3 = -^nV + ^A(v-i5../^). (A.15) 
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Therefore, the expression (jA.6[) takes the form 

+ 5P)g^, = - (l - I) (□ - aA) h 

+ yA/i. (A.16) 

Using this expression and the fact that (/'"^(JG^^ — 
\I\h/2 one finds that the sum of the first two nonlocal 
terms in (jA.14p is in fact local (here we consider all the 
contributions to (IA.14P modulo the overall factor Af^/4) 



2 □ - aA ' 



Using the relation (IA.6I) and its transpose^ one can 
show that the last three terms in (|A.14p reduce to 



J I D + P 



-25R^''^-^{S'^ + SP)g^, 



+2 {5^ + SP) -^-^ 5Gn 
J d^xg^'^ I Ah^'^SR^,, - SKhl^ + A (l - I) 

/i i .(A.18) 



1 a^A^ 
+2 a SR^" ^ - —:— h 



1 



U + P 2 U-aK 

Finally, the third term in (jA.14p can be transformed as 



-^Jd^xg^/^ g'"^ {5^ + 5P) {5^ + 5P] 



1 



d'^x g^'^ { a - 2 - 2B{A + 2B) 
2a -C 



2a 
+Ah^ " 



hOh 



C ~B{A + 2B) (8 + 2a) - - 
- B{A + 2B) (8a + 2a^ + 8^ 



1 



□ - aA 



1 



+32A^BM + 2B) h'"' — ^ h^^ 
U + P 



-4A^ Ji?""^- 



1 



U + P 



■ 5Rau 



(A.19) 



To calculate (Ji?^" x (5_R^y-terms in the above two ex- 
pressions we use another form of the Ricci tensor varia- 



® Functional transposition of 5n should take into account the vari- 
ation of the integration measure g^/^ because, as mentioned 
above, the operator g^^'^{5\I\) is not symmetric in contrast to 
symmetric g^/^D. 



tion which holds in view of (lA.lSp and ([T3 
1 



5R^iv\ 



(A)dS 



{U + P)h^,^- -Kg^^h 
-A( ^-^"j [K.-\g^uh). (A.20) 



This form contains as a whole the operator U + P which 
cancels (□ + P)^^ and renders a part of terms explicitly 
local. The result reads 



+ ^h^^g^^5G,.^ = ^ j d'xg'/^h\{k.n) "/ 



j d'^x SRt^'^ 5R^^ 
d^xg'/^ |i 

4(— l)"n^" 

4(^+l)(^-5)'''"sT7 

Using this relation and collecting in (|A.18|) and (|A.19|) the 
coefficients of two nonlocal structures we get the nonlocal 
part of (|A.14I) . In particular, the overall coefficient of the 



h^'^{U + P) ^h^^ nonlocality turns out to be 



16A2 



B{A + 2B) 



(A.22) 



because — A^ = 8i?(A + 2B). Similarly the coefficient 
of /i(n — aA)~-'^ft, equals 



-K^B{A + 2B) 4 + a + - 
' a 



B[A + 2B) 



A' C 



(A.23) 



This finally gives the nonlocal terms of Eq. p^ and ex- 
plains the origin of the effective Planck mass ([T9| . 



2. Local terms 

Local terms of the quadratic part of the action are 
contained both in (jA.13p and (jA.14p . The most eco- 
nomical way to simplify the second variations of Ricci 
and Einstein tensors in (jA.lSp is to use the identity 
(<52i?/^-)g^, - g^-{6^G,,,) = SiSR^^'^g^, - g'^'^SG^,) - 
hf^i^SR^" — h'^'^SGfj^i, and then repeat this trick for first 
order variations. In this way the second order variations 
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reduce to the calculation of 



J d^x <?i/2 5^R ^ J d^x I 2Kh% - i Kh^ 



^h^"'SR^,^- -hSR} , (A.24) 



and the first three lines of (lA.lSp simplify to 

J d^x S^^hSR+^Ah^ + SG'^'hf,^ 
-Ihl.R+^hg^JR^'' 
+ {S'R^-)g^,-gf^'^{S'G^,) 

= f d'^xQ^/^ l^h'"'Dhf,^--hDh 



-^Ahl^ + Uh^Y (A.25) 



For the calculation of the last term of (|A.13p consider 
the second order variation of the identities D^^i/ = and 

(52^)5^, = -2((53)V' (A.26) 
{S^P) g,,, = ~2dP h^, + A S^R,,, + 2B6R V 

+Bg^,6^R, (A.27) 



and again use the first order variations (jA.4[) - (IA.5l) . 
Then bearing in mind that g'"'{SD)hf,^ = S{g''''D)hf,^ + 
ht'^Dhf,^ = {6D)h - U{hl^) + h^'^Uh^^ we reduce the 
last term of S(2) to 



d'xg''^g^^{5^^ + 5'P)g,^ 

j d^x I /i D/i _ ih^'^'Uhf,^ - a S^R 
+2A(C-«) 

= i 1 d'xg'/^ { (1 - f ) ''^^ - (2 + f ) /^""^nv 

+A/.^.(2C-^)+^/.^(f-C)}. (A.28) 



Then collecting the local terms from (|A.17P - (|A.19|) 
with the SR^'^ x (5i?^y-terms given by (IA.2ip together 
with the local terms ()A.25P and (|A.28P we finally get the 
local part of (gH])- 
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